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A simple, physically motivated, scaling hypothesis, which becomes exact in important limits, 
yields estimates for the ground-state energy of large, composed, systems in terms of the ground- 
state energy of its building blocks. The concept is illustrated for the electron liquid, and the Hubbard 
model. By means of this scaling argument the energy of the one- dimensional half-filled Hubbard 
model is estimated from that of a 2-site Hubbard dimer, obtaining quantitative agreement with the 
exact one-dimensional Bethe-Ansatz solution, and the energies of the two- and three-dimensional 
half-filled Hubbard models are estimated from the one-dimensional energy, recovering exact results 
for L'^ ^ and U oo and coming close to Quantum Monte Carlo data for intermediate U. 
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Scaling techniques have propped up research since the 
early days of modern science. In the traditional version, 
which we call discrete scaling, one measures the prop- 
erties of a body of diameter D to determine those of a 
similar body of diameter D' , typically much larger than 
D. This approach is widely employed, e.g., in fluid dy- 
namics. A more recent approach, extensively explored in 
the theory of critical phenomena, focuses on the relation 
between differential increments in the size L of a complex 
system and the resulting changes in its physical proper- 
ties. To emphasize that L is a continuous variable, we 
refer to this procedure as continuous scaling. Continuous 
scaling finds a parallel in studies, also pioneered in the 
theory of critical phenomena, that treat the dimensional- 
ity d of a system as a continuous variable, and hence will 
be called continuous dimensional scaling. A classic illus- 
tration is the e expansion, which extracts the properties 
of (4 — e)-dimensional systems from the exact mean-field 
solution of the 4-dimensional Gaussian modeli Other re- 
cent examples are found in dynamical mean-field theory^ 
(DMFT) and in density-functional theory (DFT) for the 
Heisenberg modeli^ 

Short as it is, this overview points to an as yet un- 
probed form of scaling: discrete dimensional scaling. We 
are unaware of attempts to extract the properties of a 
system of dimensionality d from the known properties 
of the same system in a different dimensionality. In the 
present work we explore this possibility. As a specific 
nontrivial model, to illustrate the general idea of discrete 
dimensional scaling, we choose the Hubbard model. 

The Hubbard model was originally proposed on a 
three-dimensional lattice for the description of corre- 
lations in narrow-band materials.* Its two-dimensional 
version is often considered the correct minimal model 
for describing the cupper-oxide planes in cuprate 
superconductorSf" The one-dimensional Hubbard model 
has received much attention due to the Luttinger liquid 
and Mott insulator phases it displays i& In the limit of 



infinite number of dimensions the Hubbard model was at 
the beginning of DMFT.^ In all dimensions, the model 
serves as a laboratory and prime test case for our under- 
standing of many-body physics and, in particular, strong 
Coulomb correlations. In standard notation and for any 
dimensionality d, the model reads 

H = f+U=-tY^ clcj. + UY,4^'^^^Cy^U (1) 

where t and U are the kinetic energy and interaction 
parameters, respectively, and (ij) denotes a sum over 
nearest neighbors on a d-dimensional lattice. 

The one-dimensional homogeneous Hubbard model has 
an exact solution in terms of the Bethe AnsatZfSii one of 
the key results of which is a closed expression for the per- 
site ground-state energy at half filling, as a function of 
U, 

tN, Jo x[l + exp{Ux/2)] ' ^ ^ 

where is the number of sites and Jq and Ji are zero 
and first order Bessel functions. Such simple expressions 
for the total energy of nontrivial models are scarce in 
many-body physics, and where they exist they provide, 
in addition to considerable algebraic simplifications, in- 
sights into the nature of the model that would be hard 
to obtain from numerical data on its own. 

The purpose of the present paper is to investigate 
to which extent knowledge aquired in low-dimensional 
situations, such as Eq. ((2Jl, can be useful in higher- 
dimensional systems, where exact solutions are not avail- 
able. Mathematically, low and higher dimensional sys- 
tems are solutions of the same type of differential equa- 
tion, Schrodinger's equation. Physically, higher dimen- 
sional systems are built out of lower-dimensional build- 
ing blocks. Both of these facts suggest that higher- 
dimensional systems cannot be arbitrarily different from 
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lower-dimensional ones. Mathematical and structural 
constraints imply the existence of (perhaps complicated) 
connections between quantities pertaining to one and to 
higher-dimensional systems. Here we identify and ex- 
plore such constraints. 

A first illustration of the concept of dimensional scaling 
is provided by bound homogeneous Coulomb systems, 
such as the Fermi liquid, described by the Hamiltonian 



Hch 



w, 



(3) 



where W represents the constant potential energy of the 
uniform charge background, or of confining walls. For 
such a system satisfaction of the quantum-mechanical 
virial theoremi guarantees dynamical stability of the 
ground state. This theorem states that the total ground- 
state energy (comprising kinetic, potential and inter- 
action energy) is the negative of the kinetic energy, 
i.c.,^ Etot ~ —Ekin- Any relation of the form Etot = 
AEkirn with constant A, allows to write the identity 
Etot{d)/Etot{d') = Ekin{d)/Ekin{d'), which implies that 
the d-dimensional total energy can be obtained by scaling 
the d'-dimensional total energy according to 



Etot {d) = Ekin (d) 



Etotjd') 

Ekin{d') 



(4) 



Equation Q is a first example of dimensional scaling, 
showing that the solution of a low-dimensional problem 
can generate a solution to a higher-dimensional problem 
of the same type. The virial theorem plays the role of 
the constraint mentioned above. 

Next, we explore the question if similar expressions 
can also be obtained for the Hubbard model, whose total 
ground-state (GS) energy we denote by E{U, d). Instead 
of the virial theorem, whose above form does not apply 
to the Hubbard model, we here base our analysis on the 
Hellman-Feynman theorem, which allows to cast the ki- 
netic energy Ekin{U, t) = E{U, t) - Em{U, d) as 



Ekin{U, t) 



dt 



dU 



(5) 



The second term on the right-hand side vanishes both at 
small U, where the total ground-state energy E{U, d) is 
linear in U, and at large U, where E{U,d) oc 1/U. Hence, 
in both limits E{U, d) — +Ekin{U, d). As for the electron 
liquid, this equality gives rise to the dimensional-scaling 
expression for the d-dimensional ground-state energy. 



E{U,d) ^ E^^iU^d) = Ek^n{U,d) 



E{U, d) 

Ek,niU,d') 



(6) 



but unlike for the electron liquid, this is exact only in 
two limits, and approximate in between. Its use at all U 
embodies a linear dimensional-scaling hypothesis, whose 
validity must be carefully checked. Our first check is 
an application of Eq. (0 to one-dimensional Hubbard 
chains, whose energies we try to recover by means of 



Eq. © from that of a zero-dimensional Hubbard model. 
A truly zero-dimensional model is useless for the present 
purpose, but a two-site system is as close to zero dimen- 
sionality as one can get without losing either of the two 
defining terms of the Hubbard model, the on-site inter- 
action and the nearest-neighbour hopping. 
Applied to this case, Eq. © reads 



E{U,1) ^ E'^-^iU,!) ^ EMniU.l)' 



E{U, L = 2) 
5fc„(C/,L = 2)' 



(7) 



where L is the number of sites. i?(?7, 1) is exactly known 
from Eq. which thus allows a stringent test of the 
concept of dimensional scaling. The two-site Hubbard 
model can be diagonalized analytically for any U . For a 
half-filled band the exact per-site ground-state energy is 



E{U, L^2) 



U 



El 
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The kinetic energies in Eq. Q can be obtained from ^ 
and © via Ekm = tdE/dt. All ingredients on the right- 
hand side of lO arc thus known, and the resulting pre- 
diction for E{U, 1) can be compared to the exact result 
for the infinite chain. This comparison, in Fig. ^ shows 
that linear dimensional scaling (dotted curve) almost co- 
incides with the exact energy (dashed curve) , illustrating 
that Eq. {Tj) with © is a near-exact representation of the 
true ground-state energy, even at intermediate U. 

In more complicated cases the exact kinetic energies 
may not be known. We therefore also consider an ap- 
proximation to Eq. Q in which noninteracting kinetic 
energies are used in the scaling relation. For the dimer 
and the chain the noninteracting per-site kinetic energies 
are E{U = 0, L = 2) = and E{U = 0, d = 1) = -4t/7r, 
respectively. Hence, the noninteracting dimensional- 
scaling prediction for the per-site energy of the infinite 
chain is 



^ ' ' ^ ' 'e{0,L = 2) it 

(9) 

We note that in the limits of very weak and very strong 
interactions this equation predicts E^^{U — > 0,1) = 
-At/n and E^^{U ^ oo,l) = 0, respectively. These 
are the same values also obtained from the Bethe-Ansatz 
expression showing that in spite of the additional 
approximation entailed by using the noninteracting ki- 
netic energy ratio instead of the interacting one, the di- 
mensional scaling Ansatz is still exact for very weak and 
very strong interactions. For intermediate U , a numer- 
ical comparison between the exact Bethe-Ansatz result 
and the dimensional scaling estimate Q is presented in 
Fig. n (full curve). 

Figure n shows that, in spite of its simplicity, the di- 
mensional scaling expression ((TJ reproduces the function 
E{U,d = 1) almost exactly, while © provides a very 
reasonable approximation to it. Hence, the ground-state 
energy of the infinite chain is quantitatively recovered by 
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FIG. 1: Ground-state energy of the one-dimensional half- 
filled Hubbard model at positive U. Lower (dashed) curve: 
exact result, Eq. ||5J. Middle (dotted) curve: present expres- 
sion O , representing linear dimensional scaling with interact- 
ing kinetic energies. Upper (full) curve: present expression 
@, representing linear dimensional scaling with noninteract- 
ing kinetic energies. 



suitable scaling of a simple dimer calculation. Crucially, 
this scaling involves an energy ratio, and not a length 
ratio. 

After this exploration of the concept of dimensional 
scaling in situations in which the exact result is known 
in closed form, we now turn to a situation in which 
the higher-dimensional result is known only numerically. 
Namely, we attempt to predict the ground-state energy of 
the half-filled two-dimensional repulsive Hubbard model 
by applying virial-based dimensional scaling to the Lieb- 
Wu Bethe-Ansatz expression This expression thus 
ceases to be our target to become the input on the right- 
hand side of Eq. ©. 

The d-dimensional kinetic energy Ekin (d) of the Hub- 
bard model is already an in general unknown function 
of d, so that Eq. ^ cannot be used directly. We can, 
however, still employ noninteracting kinetic energies in 
the scaling factor. Equation © now reads 



E{U, 2) « E^'^-^'iU, 2) = E{0, 2) 



E{U,l) 
S(0,1)' 



(10) 



Note that all quantities appearing on the right-hand- 
side of (|10|l are known: E{U, 1) is given by the Lieb- 
Wu formula lO, while E{0,d) has the valuesifiJii^ 
E{0,l)/tN, = -1.2732, E{0,2)/tNs = -1.6211, and 
E{0,3)/tNs = -2.0048. As before, Eq. (dOJ, which is by 
construction exact in d = 1 for any U, and at [/ = in 
any d, is also already exact in any d in the limit U oo: 
it follows from Eq. Q that E{U ^ oo, 1) = 0, and hence 
j^LDSfjj ^ Qo cf) — 0, which is the correct U ^ oo limit 
for the half-filled Hubbard model in any d. Dimensional 
scaling with noninteracting kinetic energies is thus still 
exact both for very weak and very strong correlations. 



FIG. 2: Ground-state energy of the two-dimensional half- 
filled Hubbard model. Filled triangles: VMC data extrap- 
olated to the thermodynamic limit, from Ref. ITol Dashed 
curve: Hartree-Fock energy. Dotted curve: [/ — > cxd expan- 
sion of Ref. ITgI Full curve: present expression IIOII . 



A more comprehensive check on (|10|l is obtained 
by comparing with precise numerical data obtained by 
Monte Carlo simulations ii2*i2iiii4 Below we adopt as a 
benchmark the Variational Monte Carlo (VMC) data of 
Yokoyama and Shiba.^° Their values do not deviate sub- 
stantially from those obtained by other groupsjiS, Fig- 
ure 12 compares, for the d = 2 Hubbard model, the di- 
mensional scaling expression H10|) with these benchmark 
data, and with known large and small U limits. 

Finally, we turn to the three-dimensional case. The 3d 
energies could be constructed in terms of the 2d ones, in 
the same way as before, but since the 2d energies them- 
selves are not available in closed form this is not imme- 
diately practical. Instead, we subject the dimensional 
scaling idea to a more severe test, and try to predict the 
3d energies directly from the Id energies. Equation © 
then becomes 



E{U, 3) w E^'^-^iU, 3) = £:(0, 3) 



E{U,1) 
E{0,1)- 



(11) 



A comparison of the resulting curve with QMC data and 
known large and small U limits is shown in Fig. |31 

Figures |21 and O show how Id to 2d and Id to 3d di- 
mensional scaling, as expressed by Eqs. ((Tn|l and l(TT|l . 
merges into the exact expressions in the limits J7 — > and 
U ^ oo, and comes close to the QMC data at interme- 
diate U. Not unexpectedly, the step from d = 1 directly 
to d = 3 entails a larger margin of error than those from 
d — to d — 1 and from d — 1 to d = 2, because the di- 
mensional gap to be bridged by a simple scaling hypoth- 
esis is larger. The concatenation of the d = to d = 1 
expression with the d = 1 to d = 3 expression shows that 
an exact solution of the two-site system, followed by a 
double application of linear dimensional scaling, is suffi- 
cient to predict energies of the three-dimensional system 
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FIG. 3: Ground-state energy of the three-dimensional half- 
filled Hubbard model. Filled triangles: VMC data on a 6 x 6 x 
6 cluster (not extrapolated to the thermodynamic limit) from 
Ref. ITol Dashed curve: Hartree-Fock energy. Dotted curve: 
U ~> oo expansion of Ref. 16. Full curve: present expression 
Q, representing linear dimensional scaling. 



with a worst error of the order of 20 percent (and much 
less than that for most values of U), without empirical 
adjustments or fitting parameters. Remarkably, a two- 
site system thus contains enough nontrivial many-body 
effects that a simple scaling hypothesis is sufficient to re- 
cover from it a large part of the energy of chains, planes 
and bulk systems. 

In summary, we have introduced the concept of di- 
mensional scaling of ground-state energies. One explicit 



realization of this concept, based on the virial theorem, 
is exact for the electron hquids. Another, based on the 
Hellman-Feynman theorem, allows to predict the energies 
of an infinite Hubbard chain from that of a dimer with 
unexpected accuracy, and continues reasonable when ex- 
tended to planes and bulk systems. Although the non- 
interacting kinetic energies in the scaling factor intro- 
duce an additional approximation, the dimensional scal- 
ing ansatz remains exact in the limits of weak and strong 
interactions. 

These observations suggest the following questions for 
further investigation: (i) Does dimensional scaling also 
hold for other Hamiltonians? (ii) Is there a way to sys- 
tematically improve on the dimensional scaling expres- 
sion, reducing the remaining difference to the benchmark 
data without fitting or empirical input? (iii) Does dimen- 
sional scaling hold for other quantities than energies? A 
partial answer to question (i) is given by Eq. Q), which 
shows that dimensional scaling can also be expected for 
the electron liquids. A related, but not identical, type of 
dimensional scaling was recently also observed for the 
Heisenberg model. Preliminary work on the negative 
U (attractive) Hubbard model also indicates the exis- 
tence of simple scaling lawsiii The present approach is 
expected to be particularly useful for model Hamiltoni- 
ans that are not Bethe-Ansatz integrable, but can be 
solved in the 2-site limit, such as the periodic Anderson 
model. As for question (ii), the answer is 'yes'. A simple 
nonempirical way to improve on virial-based dimensional 
scaling will be presented in a forthcoming publication.^'^ 
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